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O \ Abstract. Recent observations by Chandra and XMM-Newton demonstrate that the central gas 

in "cooling flow" galaxy clusters has a mass cooling rate that decreases rapidly with decreasing 
' temperature. This contrasts the predictions of a steady state cooling flow model. On the basis of these 

. observational results, the gas can be in a steady state only if a steady temperature dependent heating 

mechanism is present; alternatively the gas could be in an unsteady state, i.e., heated intermittently. 
Intermittent heating can be produced by accretion on the supermassive black hole residing in the 
central cluster galaxy, via Compton heating. This mechanism can be effective provided that the 
radiation temperature of the emitted spectrum is higher than the gas temperature. Here we explore 
whether this heating mechanism can be at the origin of the enigmatic behavior of the hot gas in 
0^ ■ the central regions of "cooling flow" clusters. Although several characteristics of Compton heating 

appear attractive in this respect, we find that the fraction of absorbed heating for realistic gas and 
radiation temperatures falls short by two orders of magnitude of the required heating. 
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^ : INTRODUCTION 

X-ray observations of the central regions of a large fraction of galaxy clusters prior to 
Chandra and XMM — Newton were interpreted in terms of cooling and condensing of 
the intracluster medium (ICM), leading to a subsonic, steady central inflow called cool- 
ing flow (Cowie & Binney 1977, Fabian & Nulsen 1977). The measurements of radiative 
cooling times lower than a Hubble time within a cooling radius (r coo i) of ^ few hundreds 
of kiloparsecs were at the basis of this idea. However, cool gas in the cluster cores (in 
the form of massive-star formation or formation of low-mass stars, optical emission line 
nebulae or cold gas) has never been discovered in a quantity large enough to fit with 
the steady state cooling flow predictions (e.g., Donahue & Voit 2003). Moreover recent 
observations from the Chandra and XMM-Newton satellites have ruled out the simple 
steady state cooling flow model (e.g., Molendi & Pizzolato 2001, Peterson et al. 2003). 
In this model one expects emission from gas over the entire temperature range (from 
Tmnx, the ambient ICM temperature, down to temperature values at which the ICM is 
undetectable in the X-rays), with the same mass cooling rate M at each temperature. 
The strengths of a few observed emission lines reveal how much gas cools through 
each temperature, and high resolution spectroscopic observations now show a deficit of 
emission relative to the cooling flow predictions from gas below ~ T max /3 (Peterson 
et al. 2003). In addition, the spectra show increasingly less emission at lower tempera- 
tures than the cooling flow model would predict. Empirically, the differential luminosity 
distribution AL/AT T 1 ^ 2 , instead of being temperature-independent as expected in 



the radiative cooling flow model. Lower spectral resolution observations, as from the 
Chandra ACIS-S detector, also suggest significantly lower mass cooling rates than ob- 
tained from previous analyses of ROSAT and ASCA data (McNamara et al. 2000). 

These observational results have produced an important astrophysical puzzle and a 
number of ideas have been suggested to solve it, but none has been proven conclusive 
yet. For example, cooling may be opposed by heating (Soker et al. 2001) or thermal 
conduction may suppress cooling (Ruskowkski & Begelman 2002). 'Cooling flows' 
generally have embedded non-thermal radio sources and this association suggests that 
feedback from an AGN may help suppress the cooling of the ICM. This feedback may 
come from the impact of the radio jets on the cooling gas (Reynolds et al. 2002, Omma et 
al. 2003) or from mixing and turbulent heating as buoyant radio plasma rises through the 
ICM (Briiggen & Kaiser 2002). Another plausible (intermittent) heating source could be 
Compton heating resulting from accretion of the cooling gas on the central supermassive 
black hole of the central galaxy in the cluster core. Although the gas over the body of 
the galaxy is optically thin, numerical simulations showed that this mechanism is very 
effective in heating the interstellar medium of giant elliptical galaxies (Ciotti & Ostriker 
2001). Here we explore whether it can provide effective feedback also in the case of 
cluster cores. 

The bolometric luminosity associated with accretion in the galactic nucleus is Lbh = 
£ Mbhc 2 , where Mbh is the accretion rate on the black hole, £ is the accretion efficiency 
(usually spanning the range 0.001 ^£^0.1) and c is the speed of light. Taking Mbh of 
the order of few 10s of M yr _1 , as suggested by the most recent estimates of the mass 
accretion rate in 'cooling flows' (Peterson et al. 2003), the power available for heating of 
the ICM turns out to be ~ 1 47 erg s 1 . The power required to balance the cooling of the 
ICM in the 'cooling flow' region must be of the order of its observed X-ray luminosity; 
this goes from ~ 10 43 erg s _1 in the Virgo cluster up to ~ 10 45 erg s _1 in the most 
massive clusters like A1835. Therefore it seems that during the phases of accretion 
there could be enough power to balance the cooling. The problem is: How much of 
Lbh is actually trapped by the inflowing gas and therefore is effectively available for its 
heating? It is clear that for extremely low opacities Lbh would be unable to affect the 
flow. In the following we estimate how much of Lbh is absorbed by the ICM in the core 
regions of galaxy clusters, for typical temperature and density profiles. 



SETTING THE PROBLEM 



The number of photon-electron interactions per unit volume and in the time interval At 
at any radius r in a plasma can be written as: 



N re (v,r) = N Y (v,r) x = x ^(vKW 
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where n s (r) is the electron number density, h is the Planck constant and Okn is the 
Klein-Nishina electron scattering cross-section (Lang 1980): 
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where x = v/vj, Vj = m e c 2 /h is the Thomson frequency, m e is the electron mass and 
Or is the Thomson cross section. Here for simplicity we assume a gray absorption, i.e.: 
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^BH(v,r) =/ BH (v) xL BH (r), / /bh(v)Jv = 1, (3) 
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where Lbh(/) is the bolometric accretion luminosity that reaches the radius r from the 
nucleus. 

The gas Compton heating (or cooling) per unit frequency at radius r is given by 
AE = —Ny s (v)AEy(v,T), where AE is the internal energy per unit volume gained (or 
lost) by the gas from radiation at frequency v, and AEy is the energy variation of a 
photon of frequency v interacting with an electron of gas at temperature T. A simple 
approximation for the energy transfer factor is 1 : 
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where k B is the Boltzmann constant. After substitution of ([T]) and @ in the expression 
for the Compton heating, and integration over all frequencies, one has: 
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where nt(r) is the total number density, 
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and the spectral temperature Tq is given by: 

_ m e C 2 r x(\+X 2 )f Bll (v)OKN(v) , n , 
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where m e c 2 /4£ B = 1.48 x 10 9 K. 

Integrating on r the equation of energy conservation <9Lbh(/) /dr= —4iir 2 dE{r) / dt, 
one has 

LBH(r)=L BH (0)exp{|i^ f E {r)^f\ ( 1 - ^- ) dr) . (8) 
[3m e c z Jo n t (r) \ T(r) ) J 

From the above equation we can compute the amount of energy actually trapped 
to heat the gas, for any gas temperature and density profile and any spectral energy 
distribution of the nuclear photons. 



1 This formula reproduces the well known relations AEy ~ l.5ksT — m e c 2 x + 0(1 /x) for relativistic 
photon energy (hv » m e c 2 ), and AEy ~ Ak^Tx — m e c 2 x 2 + 0(x 3 ) in the classical limit (hv <C m e c 2 ). 



We assume that the spectral energy distribution /bh(v) is made of two distinct contri- 
butions. The first is a non thermal distribution of total luminosity Lx(v) = /x(v) X Lx, 
with 
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where hv^ is the spectrum break energy and ^ and £1 + £,2 are the spectral slopes at low 
and high frequencies. The second is a blackbody distribution Luv(v) = /uv( v ) x ^uv 
at a temperature Tyy, with 
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For these distributions fx(v)dv = 1 and Jq fuy(v)dv = 1. We also assume that 
Luv = &Lx, where M is a dimensionless parameter measuring the relative importance 
of the "UV bump" with respect to the high energy part of the spectral energy distribution. 
Following this choice, the frequency distribution of Lbh in © can be written as 



andL BH (0) =L X + L UV 
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THE RESULTS 

We first estimate the dependence of Tq and 7c, i.e., of the integrals in © and Q, 
on the free parameters Tjjv, Vb, £ji> <^2- We adopt = 0.9 and £2 — 0-7; these values 
reproduce the observed spectral shapes of the X-ray and y-ray emission of AGNs (e.g., 
Nandra & Pounds 1994; Lu & Yu 1999). The coefficient Tq is evaluated numerically by 
considering the range 0.1 < Vb/vj < 10 for fx and the classical limit for /bh(v), i.e., 
£b%v <^ hvj. We obtain 

r c = r x + ^r uv ~ ^— | x cr T . (12) 

A similar evaluation of the dimensionless integral on the r.h.s. of Q gives 

^ 8.6xl0 7 x(Vb/vT) ai96 + ^r uv 
C " 0.77x(v b /v T )- ao6 + ^ 1 j ' U; 

For observed values of v b /v T ~ 0.2 and @ ~ 1 (e.g., Fabian 1996), T c = 4 x 10 7 K 
(independently of ?uv)- This value for 7c is very similar to that derived by Sazonov et 
al. (2003) from composite mean QSOs spectra. 

We next performed the integration in © to estimate Lbh(^) — ^bh(O), which is the 
amount of energy actually trapped by the ICM within a radius r. The integration was 
extended out to the 'cooling radius' (r coo iX within which a 'cooling flow' could develop. 




FIGURE 1. The amount of nuclear radiation that actually heats the ICM within r coo ]. Vb/vx is fixed at 
0.2. The upper set of curves refers to Tq — T = 5 x 10 K, the lower set to 7c — T = 10 7 K. Solid lines 
refer to M = 1, dashed ones to M = 2 and dotted ones to S% = 0.5. 



For low 7c values (~ 1 — 2 keV, as allowed by the study of Sazonov et al. 2003), 
r coo \ encloses the only region within which the ICM temperature T is lower than Tq. 
We assume that the gas is isothermal at a temperature T within r coo i. The adopted 
ICM density profile is a deprojection of the /3— model commonly used to reproduce 
the observed X-ray surface brightness (i.e., Sarazin 1986): p gas = po[l + (r/r c ) 2 ]~ 3 ^/ 2 , 
where r c is the core radius (~ 250 kpc on average) and po is the total particle density at 
r = 0. The integration in ® gives: 
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where hq = po//im p with fim p the average particle mass. For our choice of r coo i<r c , the 
function / w r coo \/r c and is independent of p\ This approximate value for the integral is 
within 25% of the true value, for /3 = 0.3 — 0.9. In Fig. 1 we plot the amount of nuclear 



radiation that actually heats the ICM within r coo i as a function of riQr coo i. The range of 
variation for nor coo i is chosen by considering the observed ranges for r coo i (50 — 200 
kpc, Peterson et al. 2003) and hq (2 x 10~ 3 —0.1 cm -3 ). Taking n e /n t = 0.5, the free 
parameters left in the expression for Lbh(/cooi) m CU are V b/ V T> & an d Tq — T . In 
Fig. 1 we consider a possible range for M = 0.5 — 2 and two extreme cases for the 
difference Tq — T, corresponding to a high Tq value [~ (5 — 6) x 10 7 K] and a low one 
(I c ~2x 10 7 K). Note that T c varies by < 10% for v b /v T = 0. 1 - 1, and so the results 
are totally unaffected by variations of Vb/vx in this range. 



CONCLUSIONS 

From Fig. 1 it appears that the power actually available for heating of the ICM within 
the cooling region goes from fewx 10~ 6 Lbh(0), in the small clusters, up to 10~ 4 Lbh(0) 
in the most massive ones. This makes Compton heating fall short by ~two orders of 
magnitude of the required heating and therefore an unplausible mechanism to balance 
the cooling of the gas in the cluster core, at variance with the situation in elliptical 
galaxies (Ciotti & Ostriker 2001). 
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